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Secondary  instability  of  a  temporally 
growing  mixing  layer 

by  Ralph  W.  Metcalfe,  Steven  A.  Orszag,  Marc  E.  Brachet,  Suresh  Menon, 

and  James  J.  Riley 

Flow  Research  Company,  Kent,  Washington  98032 

The  three-dimensional  stability  of  two-dimensional  vortical  states  of  planar 
mixing  layers  is  studied  by  direct  numerical  integration  of  the  Navier-Stokes 
equations.  Small-scale  instabilities  are  shown  to  exist  for  spanwise  scales  at 
which  classical  linear  modes  are  stable.  These  modes  grow  on  convective  time 
scales,  extract  their  energy  from  the  mean  flow,  and  persist  to  moderately  low 
Reynolds  numbers.  Their  growth  rates  are  comparable  to  the  most  rapidly  growing 
inviscid  instability  and  to  the  growth  rates  of  two-dimensional  subharmonic 
(pairing)  modes.  At  high  amplitudes,  they  can  evolve  into  pairs  of  counter¬ 
rotating,  streamwise  vortices,  or  "ribs",  in  the  braids,  which  are  very  similar 
to  the  structures  observed  in  laboratory  experiments.  The  three-dimensional 
modes  do  not  appear  to  saturate  in  quasi-steady  states  as  do  the  purely 
two-dimensional  fundamental  and  subharmonic  modes  in  the  absence  of  pairing. 

The  subsequent  evolution  of  the  flow  depends  on  the  relative  amplitudes  of  the 
pairing  modes.  Persistent  pairings  can  inhibit  three-dimensional  instability 
and,  hence,  Keep  the  flow  predominantly  two-dimensional.  Conversely,  suppres¬ 
sion  of  the  pairing  process  can  drive  the  three-dimensional  modes  to  more 
chaotic,  turbulent-like  states.  An  analysis  of  high-resolution  simulations  of 
fully  turbulent  mixing  layers  confirms  the  existence  of  riblike  structures  and 
that  their  coherence  depends  strongly  on  the  presence  of  the  two-dimensional 
pairing  mode:; . 


1 .  Introduce  ion 

Free  shear  flows,  like  those  of  mixing  layers  and  jets,  differ  from  wall- 
bounded  flows  in  that  they  typically  have  inflexional  mean  velocity  profiles 
and,  hence,  are  subject  Co  inviscid  instabilities.  Thus,  it  may  be  thought 
that  the  process  of  transition  to  turbulence  in  free  shear  flows  would  be 
directly  amenable  to  analysis.  Indeed,  observations  by  Winant  &  Browand 
(1974);  Brown  &  Roshko  (1974);  Wygnanski,  Oster,  Fiedler  &  Dziomba  (1979); 

Ho  &  Huang  (1982);  Hussain  (1983a),  and  others  show  the  central  role  played  by 
two-dimensional  dynamic  processes,  at  least  through  transitional  regimes,  in 
these  flows.  While  three-dimensional  small  scales  are  observed  (Mik3ad  1972; 
Bernal,  Breidenthal,  Brown,  Konrad  &  Roshko  1979),  they  may  not  necessarily 
destroy  the  large-scale  two-dimensional  structures  (Browand  &  Troutt  1980).  In 
contrast,  studies  of  wall-bounded  flows  have  emphasized  the  central  role  of 
three-dimensional  effects  in  the  breakdown  to  turbulence. 

In  this  paper,  we  investigate  the  interaction  between  linear  and  nonlinear 

two-  and  three-dimensional  flow  states  that  can  arise  during  the  early  stages 

of  evolution  of  a  temporally  growing  turbulent  mixing  layer.  It  is  shown  that 

certain  two-dimensional,  nonlinear  states  (coherent,  spanwise  vortical  modes) 

are  strongly  unstable  to  small,  three-dimensional  perturbations,  and  that  these 

perturbations  can  evolve  into  strearawise,  counterrotating  vortices  similar  to 

those  observed  experimentally  (Bernal  1981)  anc  modeled  analytically 

(Pierrehumbert  4  Widnall  1982;  Lin  &  Corcos  1934).  We  find  that  the  two-  or 

« 

three-dimensional  character  of  the  mixing  layer  depends  crucially  on  the 
initial  conditions,  as  there  is  a  close  competition  between  the  various  modes 
of  instability. 

The  approach  followed  here  is  similar  to  that  used  by  Orszag  &  Patera  (198G, 
1983)  in  their  study  of  secondary  instabilities  in  wall-bounded  flows.  The 


parallel  laminar  flow  is  perturbed  initially  by  either  a  linear  or  a  finite- 
amplitude  two-dimensional  disturbance  that  is  allowed  to  evolve  and  to  saturate 
in  a  quasi-steady  state.  The  stability  of  this  finite-amplitude  vortical  state 
to  both  subharmonic  (pairing)  two-dimensional  modes  and  smaller-scale  three- 
dimensional  modes  is  then  studied  by  numerical  solution  of  the  full  three- 
dimensional  time-dependent  Navier-Stokes  equations.  To  relate  these  simulations 
to  the  evolution  of  a  turbulent  mixing  layer,  we  also  examine  the  interaction 
between  the  evolving  two-dimensional  modes  in  their  linear  and  nonlinear  states 
and  a  broad-band,  three-dimensional  background  noise  spectrum. 

The  character  of  the  pairing  instability  was  first  explained  theoretically 
by  Kelly  (1967)  and  numerically  by  Patnaik,  Sherman  &  Corcos  (1976)  and  Collins 
(1982)  for  stratified  flows,  and  by  Riley  &  Metcalfe  (1980)  and  Pierrehumbert  & 
Widnall  (1932)  for  unstratified  flows.  The  nature  of  the  two-dimensional 
vortical  pairing  as  well  as  a  model  for  streamwise  vortical  motion  have  been 
investigated  numerically  and  theoretically  (Corcos  &  Sherman  1984;  Corcos  &  Lin 
1984;  Lin  &  Corcos  1984).  Experimentally,  coherent  pairing  of  large-scale 
vortical  structures  in  turbulent  mixing  layers  at  high  Reynolds  numbers  was 
identified  by  Brown  &  Roshko  (1974).  Significant  secondary  three-dimensional 
instabilities  in  these  flows  have  been  observed  by  Breidenthal  (1981)  and  Bernal 
(1981).  The  importance  of  these  instabilities  and  their  sensitivity  to  up¬ 
stream  perturbations  has  been  demonstrated  experimentally  by  Hussain  &  Zanan 
(1978),  Oster  &  Wygnanski  (1982),  and  Ho  &  Huang  (1982)  among  others.  There  is 
an  excellent  and  comprehensive  feview  of  the  very  extensive  literature  on  this 
topic  by  Ho  4  Huerre  (1984). 

Pierrehumbert  &  Widnall  (1982)  examined  the  linear  two-  and  three- 
dimensional  instabilities  of  a  spatially  periodic  inviseid  shear  layer  in  a 
study  closely  related  to  the  present  one.  Thev  considered  the  stability 


characteristics  of  the  model  family  of  two-dimensional  vortex-modified  mixing 
layers  with  velocity  fields, 


u  ■  sinh  z/[cosh  z  -  p  cos  x] 
w  ■  -  p  sin  x/[eosh  z  -  p  cos  x] 


(1.1) 


(Stuart  1967)  for  0  <  1,+  and  studied  subharmonic  pairing  instabilities  and 

a  "translative"  three-dimensional  instability.  In  contrast,  we  consider  here 
both  the  linear  and  nonlinear  stability  characteristics  of  time-developing 
viscous  shear  layers.  The  three-dimensional  secondary  instability  we  study  is 
both  the  analog  of  the  translative  instability  and  a  generalization  of  the 
instability  analyzed  by  Orszag  &  Patera  (1983)  for  wall-bounded  flows.  In  the 
nonlinear  state,  this  instability  is  manifest  as  the  streamwise ,  counterrotating 
vortices  (Bernal  1981)  or  "ribs"  (Hussain  1983a)  seen  in  laboratory  experiments. 


2.  Numerical  methods 

The  Navier-Stokes  equations  are  solved  in  the  form 

3v  "*  _  _2 

■jj'  ■  V  X  111  -  Vlt  ♦  .V  V  V 

V  •  v  ■  0  , 


(2.1) 

(2.2) 


t  Note  that  for  p  «.^  the  basic  flow  state  (1.1)  is  of  the  form 

tanh  z  x  ♦  p  Re  [e1  V(z)].  At  wavenumber  a  *  1,  there  are  no  fundamental 

two-dimensional  instabilities  that  can  compete  with  the  subharmonic  (  a  =  1/2) 
and  secondary  instabilities.  This  flow  state  is  an  inviscid  neutrally 
stable  perturbation  of  the  mixing  layer  tanh  z  x.  In  contrast,  the  results 
to  be  reported  in  §3  involve  unstable  fundamental  perturbations  to  the 
mixing  layer. 
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where  a)  *  V  x  v  is  the  vorticity  and  ti  =  p  ♦  1  2  ivi  is  the  pressure  head. 
Periodic  boundary  conditions  are  applied  in  the  streanwise,  x,  and  spanwise,  y, 
directions,  where 

v(x  +  47T/a,y,z,t)*v(x,y,z,t)  (2.3) 

v(x,y  ♦  2  tt/  Q ,  z  ,  t )  *  v(x,y,z,t)  , 

A 

while  the  flow  is  assumed  quiescent  (v  U^x;  U+  constants)  as  z  •+  ±“.  Note 
that  the  assumed  periodicity  length  B  is  4tr/a  (or  87 r/a)  to  accommodate  both  the 
fundamental  mode,  with  x-wavenumber  Ct,  and  its  subharmonic,  with  x-vavenumber 
a/2  (or  (V4),t 

Our  simulations  are  of  a  temporally  growing  mixing  layer.  By  avoiding  the 
requirement  of  imposing  inflow-outflow  boundary  conditions,  which  is  essential 
in  simulations  of  a  spatially  growing  flow,  a  faster,  more  efficient  code  can 
be  written.  Thus,  the  temporally  evolving  mixing  layer  can  be  simulated  at 
higher  Reynolds  numbers  and  with  better  resolution  than  the  spatial  flow  for  a 
given  level  of  computer  resources.  As  will  be  shown  later,  there  are  very 
important  linear  and  nonlinear  dynamic  features  that  are  common  to  the  two 


+  Pierrehumbert  &  Widnall  (1982)  point  out  that  Floquet  theory  implies  that  the 
Navier-Stokes  equations  linearized  about  jj_  flow  periodic  in  x  admit  solutions 
of  the  more  general  form  v(x,y,z)  ■  e1-**  V(x,y,z),  where  V  is  periodic  in  x 
with  the  same  periodicity  as  the  basic  flow  and  y  is  arbitrary.  However, 
Pierrehumbert  &  Widnall  consider  only  the  subharmonic  and  fundamental  cases. 
The  analysis,  which  has  not  yet  been  done  for  more  general  y,  may  be  able 
to  address  more  precisely  such  phenomena  as  the  "collective  interaction" 
described  in  experiments  by  Ho  &  Nosseir  (1981).  Indeed,  Busse  &  Clever 
(1979)  point  out  the  importance  these  general  y-modes  in  Benard  convec¬ 
tion.  The,  present  study  is  restricted  to  y  being  a  half-integer  multiple 
of  the  fundamental  wavenumber,  because  our  code  is  fully  nonlinear  with  t!v* 
periodicity  condition  (2.3).  Numerical  simulations  (Corcos  4  Sherman  lq8^.) 
with  values  of  y  different  from  ours  indicate  that  the  longest  wave  allowed 
by  the  grid  will  eventually  dominate,  although  the  details  of  the  "pairing" 
process  may  differ. 


WT 


flows,  so  a  detailed  analysis  of  numerical  simulations  of  a  temporal ly  growing 
flow  can  yield  important  insight  into  the  evolution  in  the  spatial  case.  There 
are  also  significant  differences  between  the  flows,  which  should  be  addressed 
by  future  simulations  using  inflow-outflow  boundary  conditions. 

We  have  used  two  independently  written  numerical  codes  for  the  simulations 
described  in  this  paper.  In  the  first,  the  dynamical  equations  are  solved  using 
pseudospectral  methods  in  which  the  flow  variables  are  expanded  in  the  series 


v(x,y,z,t) 


£ 

...1, 


u(m,n,p, t) 


inset  in/3y  _  ,  _ . 
e  e  7  T  (Z), 
P 


(2.4) 


where  n  and  p  are  integers,  while  m  is  a  half-integer  when  one  pairing  is 
allowed  and  a  whole  integer  if  all  pairing  modes  are  excluded.  Here  Z  *  f ( 2 ) 
is  a  transformed  z-coordinate  satisfying  Z  *  ♦  1  when  z  ■>  ®.  Two  choices 

of  f(z)  have  been  studied,  viz. 


Z  *  tanh  rr  (  fz!  <  ®, 
L 

[Zl 

<  1) 

(2.5) 

Z  -  - - -  (  ;  z  I 

<  ®. 

IZl  <  1), 

(2.6) 

r 5  n 

V2 


where  L  is  a  suitable  scale  factor.  With  these  mappings,  derivatives  with 
respect  to  z  are  evaluated  pseudospectraily  using  the  relations 


l 


for  (2.5)  and  ( 2.6)  , 
Time  stepping  is 
terms  are  marched  in 


i  <;  -  z:>  & 

oz  u  a.Z 


3v  1  1  2  9v 

2z  "■  I  *  ‘  Z  3Z 

respectively. 

done  by  a  fractional  step  method  in  which 
time  using  a  second-order  Adams-Bash forth 
'.-.posed  implicitly  using 


(2.7) 

(2.8) 


the  nonlinear 
scheme  while 
Crank-N icol son 


pressure  head  and  viscous  effects  are 


d  i  t  f  erenc  ing  .  This  scheme  is  globally  sec  and -or  Jet — accurate  in  t  ime ,  despite 


time  splitting  (Orszag,  Deville  &  Israeli  1985),  because  the  various  split 
operators  commute  in  the  case  of  quiescent  boundary  conditions  at  :  5  ♦  ®. 

There  is  one  further  technical  detail  regarding  this  numerical  method  that 
should  be  discussed  here.  Various  Poisson  equations,  like 

2 

—  -  (m2  ♦  n2>  n  «  g(z)  (Is!  <  «)  ,  (2. 

dzz 


.  .  n  2 

are  solved  by  expansion  in  the  eigenfunctions  of  d“/dz  : 


dz 


-  e^(z)  *  X^e^tz)  (Izl  < 


Thus,  if 


then 


IKz) 


k=*0 


V(" 


2  2?  VJ> 

♦n  ) 


(2.1! 


(2.1 


We  remark  that  this  technique  gives  spectrally  accurate  solutions,  despite  the 
fact  that  the  continuous  version  of  the  eigenvalue  problem  (2.10)  has  only  a 
continuous,  and  hence  singular,  spectrum.  Also,  note  that  all  the  eigenvalues 
X,^  are  real  and  nonpositive;  for  both  mappings  (2.5)  and  (2.6),  there  are 
precisely  three  zero  eigenvalues  X^,  X^ ,  X^.  One  of  these  zero  eigenmodes 
is  physical,  viz.  e.(z)  =  1,  but  tha  other  two  are  highly  oscillatory  and  un- 

L  • 

physical.  Indeed,  since  the  spectral  (Chebvshev)  derivative  of  T  (2)  vanishes 

P 

•»  i 

except  at  2  *  ♦  1 ,  e  (2)  *  T  (2)  is  a  zero  e i genf unc t ion  of  d“/dz“;  T  (2.)  **  (■ 
-  2  p  p  J 

at  the  Chebvshev  collocation  points  2.  *  cos  Tij/P.  The  third  zero  eigenraode 
oscillates  and  grows  roughly  like  z .  When  m  -  n  *  0,  the  incon-.pres  s  ib  i  1  i  t  v 


-8- 


constrainc  (2.2)  requires  that  this  mode  of  the  z-velocity  field  vanish  identi¬ 
cally  so  there  is  no  difficulty  with  the  zero-pressure  eigenvalues  , A, , . 

Comparisons  of  the  behavior  of  linear  Orr-Sotrrmerfeld  eigennodes  obtained 
using  mappings  (2.5)  and  (2.6)  show  that  (2.6)  gives  a  superior  representation 
of  these  modes  unless  L  is  fine  tuned,  which  is  not  convenient  in  the  nonlinear 
dynamic  runs.t  Some  representative  results  are  given  in  Table  1.  Notice 
that  as  a  increases,  the  optimal  choice  of  map  scale  L,  decreases.  Also, 
notice  that  the  accuracy  of  the  eigenvalue  is  much  more  sensitive  to  L  for  the 
hyperbolic  tangent  mapping  (2.5)  than  for  (2.6). 

This  nonlinear,  time-dependent  Navier-Stokes  code  has  been  tested  for  the 
generalized  Taylor-Green  vortex  flow  (2.12)  and  also  for  linearized  eigen¬ 
function  behavior,  with  satisfactory  agreement  being  achieved  with  power  series 
in  t  (Brachet  et  al.  1983)  and  linear  behavior,  respectively. 

The  second  code  used  in  the  simulations  is  similar  to  the  one  just 
described,  except  that  sine  and  cosine  expansions  in  z  were  used  instead  of 
Chebyshev  polynomials  (2.i).  Thus,  the  transverse  domain  extent  is  finite, 
and  care  must  be  taken  to  identify  possible  interference  effects.  Comparisons 
with  results  from  the  first  simulations  as  well  as  simulations  performed  on 
varying  size  domains  has  verified  that  such  effects  are  small  for  the  cases 


+  There  is  one  case  in  which  it  seems  that  the  hyperbolic  tangent  mapping  (2.5) 
is  more  convenient  than  the  algebraic  mapping  (2.6).  This  flow  is  the  gen¬ 
eralized  Taylor-Green  vortex  flow  that  develops  from  the  initial  conditions 

•> 

u(x,y,z,C)  =  sin  x  cos  y/cosh"  z 

.  "h 

v ( x , y , z  , 0 )  =  -  cos  x  sin  y./cosh“  z  (2.12) 

w(x,v,z,0)  =  0 

The  evolution  of  this  flow  seems  best  studied,  either  by  power  series  or  ini¬ 
tial  value  methods,  using  (2.5)  with  L  *  1.  The  time  evolution  of  this  free 
shear  flow  is  remarkably  similar  to  that  of  the  periodic  Tavlor-Creen  vortes 
(Brachet,  Meiron,  Orszag,  Nickel,  Morf  f.  Frisch  1981). 
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presented  here, 
the  size  of  the 
similar  to  that 


Aside  from  details  of  the 
arrays  used  in  the  coraputa 
described  in  Orszag  &  Rao 


input,  output 
t ions ,  the  ba 


(1974). 


buffering  routines  md 
ic  code  structure  is 


3.  Two-dimensional  instabilities 


In  this  and  the  following  sections,  results  are  reported  for  the  evolution 
of  initial  velocity  fields  of  the  form 


A10V10(2)e 


ia(x+9) 


♦  -  r 

v(x,y  ,z,0)  *  Ug(z)  x  +  Re 

^  \  i<S(x+$)  +  il3y 

*  ‘ll’ll  *  ' 


A.  v.  (z)e 


yiOx 


(3.1) 


where  0  is  the  phase  shift  between  the  fundamental  and  subharmonic  modes,  and 

<j>  is  the  corresponding  phase  shift  of  the  spanwise  mode.  The  laminar  mean 

profile  is  assumed  to  be  Uq(z)  ■  Uq  tanh  z/6^,  an  approximation  to  the  mixing 

layer  profile,  and  v..(z)  is  normalized  so  that  max  Iv..(z)l  *  1.  Here,  6.  is 
1  r  ij  Z  1J  i 

the  initial  mean  vorticity  thickness. 

The  initial  functions  v^(z)  are  normally  chosen  as  the  most  unstable 

eigenfunctions  of  the  linear  Orr-Somaer feld  equation  for  the  appropriate  wave- 

numbers  given  in  (3.1)  (Michalke  1964). t  In  this  representation,  A  is 

the  amplitude  of  the  fundamental  two-dimensional  component,  A  is  the 

1  /  ^ ,  u 

amplitude  of  its  subharmonic  or  pairing  mode,  and  A^  is  the  amplitude  of  the 
primary  three-dimensional  wave  with  a  spanwise  wavelength  equal  to  that  of  the 


+The  Reynolds  numbers  of  the  flows  discussed  below,  while  modest,  are  much 
greater  than  that  of  the  onset  of  linear  instability  (Rcr{t  =  4) ,  so  that 
even  the  linear  modes  are  effectively  inviscid.  In  this  case,  damped  modes 
may  lie  only  in  the  continuous  spectrum  (Drazin  &  Reid  1981)  and  so  are 
singular.  Whenever  equation  (3.1)  calls  for  such  a  singular  contribution  to 
the  initial  condition  (3.1),  we  choose  instead  the  flow  component  vnr  7  v;o 
of  the  fundamental  mode  (with  u,m  and  vnfa  determined  by  incompressibility) . 


L-  V. 
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furuiaraental  two-d inens ion  1 1  mode.  Time  is  n  nd imens iona 1 i zed  by  U^/ 6.  and 
space  scales  by  6^.  The  initial  co..u;tions  are  typically  chosen  so  that 

^1/2  O’  ^11  K<  Al0'  and  A10  *  0.25.  The  Reynolds  number  for  the  undisturbed 
flow  is  R  «  Un5./v. 

0  l 

In  the  absence  of  subharmonic  and  three-dimensional  perturbations 
(A1/2  g  *  A^  ■  0),  the  two-diraensionally  perturbed  flow  quickly  saturates  to 
quasi-steady  state.  In  figure  1,  a  plot  is  given  of  the  time  evolution  of  the 
two-dimensional  disturbance  energy  E^(t)  for  various  initial  amplitudes 
A^q.  The  value  for  a  is  taken  to  be  0.4446,  which  is  the  wavenumber 
corresponding  to  the  largest  growth  rate  predicted  by  linear  theory  (Michalke 
1964).  [The  range  of  inviscidly  unstable  wavenumbers  for  the  tanh  z  profile 
is  0  <  a  <  1 . ]  Here 


7 


E  <t)  “  /  lvM{**t)l  dz  • 

mn  /  mn 


where 


r  .*• 

v  ( 2  , t )  ■  >  u(n,n,p,t)  T  (z) 

mn  p 

p=*0 


(3. 


(3. 


and  u  is  defined  by  (2.4).  It  is  apparent  that  E.n  saturates  into  a  finite- 
amplitude  vortical  state  on  a  time  scale  of  order  10.  The  independence  of  the 
peak  saturation  amplitude  to  the  initial  excitation  amplitude  except  for  every 
high  initial 'amplitudes ,  which  is  evident  in  figure  1,  has  also  been  observed 
experimentally  by  Freymuth  (1966).  If  we  define  the  growth  rate  J  as 

C  *  (dE/dt ) / ( 2E)  ,  (3.. 

then  for  linear  disturbances  of  the  form 


'V  -  Re[<p(z)eia(jc_=c)] 


(  :. 
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we  have  0  *  Qc^.  The  peak  growth  rate  of  the  most  unstable  linear  mode 
( Cl  =  0.4446)  is  0  =  0.19.  Our  growth  rates,  which  are  given  later  in  the  text, 
should  be  divided  by  a  factor  of  2  for  comparison  with  the  linear  results  for 
ac^  given  by  Michalke  (1964)  to  reflect  the  different  choice  of  Uq. 

Figure  2  is  a  plot  of  the  instantaneous  spanwise  vorticity  distribution  in 
the  developed  two-dimensional  flow  for  one  of  the  runs  shown  in  figure  1.  Note 
that  while  rollup  is  occurring  by  t  *  8,  in  the  absence  of  a  subharmonic  mode, 
A^/2  q,  pairing  does  not  take  place,  and  the  flow  evolves  into  the  nonlinear, 
quasi-equilibrium  state  shown  in  figure  2b.  This  absence  of  pairing  is  analo¬ 
gous  to  that  artificially  induced  by  upstream  forcing  in  experiments  by  Miksad 
(1972),  Hussain  &  Zaraan  (1978),  Oster  &  Wygnanski  (1982),  Ho  &  Huang  (1982)  and 
others.  In  these  experiments,  the  forced  mode  is  amplified  without  also  ampli¬ 
fying  its  subharmonic.  Thus,  rollup  of  the  forced  mode  is  achieved  without 
pairing,  creating  a  region  in  the  flow  characterized  by  large-scale,  spanwise- 
coherent,  nonpairing  modes.  This  produces  countergradient  momentum  fluxes, 
interruption  in  the  growth  of  the  mixing  layer  thickness,  and  a  reversal  in  the 
sign  of  the  Reynolds  stresses  (Riley  &  Metcalfe  1980).  That  these  phenomena 
are  also  observed  experimentally  makes  this  nonlinear,  quasi-equilibrium  state 
of  interest  in  analyzing  the  physics  of  the  laboratory  flows. 

The  saturated  two-dimensional  flow  state  discussed  above  can  be  unstable  to 
subharmonic  perturbations,  A^,,  ^  in  (3.1),  for  suitable  3  (Kelly  1967).  In 
figure  3,  we  plot  the  evolution  of  the  subharmonic  perturbation  energies 

E g(t)  as  well  as  the  fundamental  two-dimensional  energy  E^(t).  Here  we 

-4 

choose  Ajg  *  0.25  and  q  *  3x10  .  In  figure  3,  the  phase  difference  9 

between  the  two  modes  is  tt/2q,  so  that  pairing  occurs.  Values  of  9  other  than 
integral  multiples  of  tt  result  in  pairing,  while  pairing  is  temporarily  inhab¬ 
ited  when  9  is  near  NTT,  with  N  an  integer.  The  cases  9  =  Nt  ar-*  an nr.alous. 


resulting  in  the  "shredding  interaction"  (P.itnn.*.  et  a  1  .  197>,/  which  is  rarely 
seen  experiraentally.t  This  subharmonic  instability  of  the  saturated  two- 
dimensional  vortical  states  is  inviscid  in  character,  as  its  growth  rate  asymp¬ 
totically  approaches  a  finite  limit  as  R  increases.  The  growth  rate  Q 

of  Che  amplitude  of  the  subharmonic  mode  is  quite  significant;  at  R  *  200, 

°l/2  0  S  when  a  “  0*A  while  q  s  0.2  for  Cx  3  0.8.  These  are  not  signi¬ 

ficantly  different  from  the  corresponding  linear  inviscid  growth  rates  of 
Orr-Sommerfeld  modes  (Miehalke  1964),  which  are  q  =  0.14  and  0.19, 

respectively. 

The  evolution  of  the  spanwise  vorticity  distribution  during  pairing  insta¬ 
bility  is  shown  in  the  contour  plots  in  figure  4  (from  Riley  &  Metcalfe  1980). 
There  is  strong  similarity  between  these  figures  and  flow  visualizations  of 
mixing  layers,  such  as  chose  of  Winant  &  Browand  (1974).  Note  that  in  two 
spatial  dimensions,  the  lines  of  constant  vorticity  act  as  fluid  markers  much 
like  the  dye  used  in  experiments.  The  evolution  of  the  pairing  process  is 
strongly  dependent  on  the  relative  initial  amplitudes  of  the  unstable  modes. 
When  q  =  A10’  the  ^undansental  node  rolls  up  first  due  to  its  higher 

growth  rate  and  shorter  saturation  time  scale.  The  subharmonic  continues  grow¬ 
ing  after  the  fundamental  saturates,  during  which  time  the  vortex  cores  gener¬ 
ated  by  the  rollup  of  the  fundamental  are  merged  into  the  subharmonic  core.  In 
this  simulation  (figure  4),  the  subharmor.ic  mode  ^  will  become  saturated 

after  about  t  ■  24,  since  there  is  no  second  subharmonic  mode  (A.,,  n)  with 

l  /  ■+  ,  u 

which  it  can  pair.  • 

t  See  Riley  4  Metcalfe  (1930)  or  Ho  4  Huerre  (1934,  p.  382)  for  additional 
plots  and  a  more  detailed  discussion  of  pnase  differences. 
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Some  insight  into  the  dynamics  oi  this  process  can  he  gamed  Iron  comparing 
figures  4  and  5.  Figure  5  is  a  series  of  velocity  vector  plots  at  the  same 
times  as  those  in  figure  4.  At  t  ■  8,  the  vortical  motion  of  the  fundamental 
cores  dominates  the  flow  and  there  are  four  nodal  points:  two  at  the  centers 
of  the  vortex  cores  and  two  stagnation  points  at  the  centers  of  the  braids.  By 
t  »  24,  the  vorticity  peaks  corresponding  to  the  fundamental  are  still  present 
(figure  4c),  although  they  are  no  longer  as  dynamically  significant  as  earlier 
(cf.  figures  5a  and  5c).  At  this  point,  there  are  basically  two  nodal  points 
in  the  flow:  one  at  the  center  of  the  subharmonic  core  and  the  other  in  the 
region  of  high  strain  at  the  center  of  the  braid.  By  t  *  32,  the  subharmonic 
mode  has  saturated  and  the  vortex  core  collapsed  into  a  quasi-equilibrium, 
nonlinear  state  like  that  for  the  unpaired  fundamental  (figure  2b). 

The  importance  of  the  pairing  process  to  the  dynamic  evolution  of  the  flow 
is  clearly  demonstrated  by  examining  plots  of  the  Reynolds  stresses  (figure  6). 
From  the  dynamic  equation  for  the  integrated  mean  kinetic  energy  E^, 


411/3 

d£M  _  [  3U 

dt  J  3z 


uw  dz 


(neglecting  viscous  effects),  it  can  be  seen  that  a  positive  Reynolds  stress 
implies  gradient  momentum  transport  (since  r’J/3z  >  0)  and  a  feed  of  energy 
from  the  mean  flow  into  the  perturbation  field,  while  negative  Reynolds  stress 
produces  countergradient  momentum  flux  and  feeds  energy  from  the  fluctuations 
back  into  the  mean.  During  the  early  stages  of  pairing,  the  Reynolds  stress  is 
predominantly  positive.  However,  by  t  =  24  (figure  be),  large  negative  regions 
(denoted  by  dashed  lines  in  the  figures)  have  developed,  and  by  t  *  32 
(figure  bd),  the  Revnolds  stress  has  become  predominantly  negative.  The 


relevance  of  the  pairing  mode  is  summarized  in  the  results  presented  in 


figure  be.  Without  the  subharmonic  pairing  mode,  the  Reynolds  stress  changes 
sign  by  t  =  16.  With  the  subharmonic  present,  however,  the  net  Reynolds  stress 
is  still  positive  at  t  *  16.  Thus,  the  suppression  of  pairing,  whether  caused 
by  forcing,  as  in  laboratory  experiments  (Hussain  6  Zaraan  1978;  Oster  & 
Wygnanski  1982;  Ho  &  Huang  1982),  or  by  eliminating  the  subharmonic  mode,  as  in 
the  numerical  simulations,  is  associated  with  a  reversal  in  the  sign  of  the 
Reynolds  stress. 

This  phenomenon  is  apparent  in  figure  3,  in  which  the  energy  in  the  funda¬ 
mental,  E10,  reaches  saturation  by  about  t  ■  10.  At  this  point,  the  Reynolds 
stress  changes  sign,  and  a  countergradient  momentum  flux  develops.  In  the 
absence  of  other  disturbances,  this  flow  will  then  evolve  into  an  oscillatory 
state  characterized  by  an  alternating  energy  exchange  between  the  mean  flow  and 
the  perturbation  fields. 

In  terms  of  turbulence  models,  the  presence  of  pairing  is  essential  to 
maintain  the  positivity  of  transport  coefficients  (such  as  eddy  viscosity). 
Since  the  eddy  viscosity,  Veddv.  is  related  to  the  Reynolds  stress  by 


3U 

Sz. 


(3.7) 


it  follows  that  suppression  of  the  pairing  corresponds  to  a  negative  eddy 
viscosity.  This  suggests  that  accurate  simulations  of  flows  with  inhibited 
pairing  may  require  the  direct  calculation  of  large-scale  structures. 

The  energetics  of  the  pairing  instability  is  revealing.  Energy  transfers 

% 

to  and  from  the  subharmonic  mode  may  be  decomposed  as 


1 


V  * 

M 


Y 


(3.3) 


where  v  involves  the  nonlinear  interaction  of  the  subharnotiic  mode  with  the 
M 

mean  flow,  the  nonlinear  interaction  of  the  subharmonic  mode  and  all 

other  two-dimensional  modes,  and  Yy  the  viscous  dissipation  of  pairing  energy. 

Here  v  and  v_  involve  sums  over  nonlinear  terms  in  the  Navier-Stok.es 
'M  ‘2D 

equations  but  are  unaffected  by  pressure;  is  proportional  to  the  enstrophy 
in  the  subharmonic  mode.  In  figure  7,  a  plot  is  given  of  these  transfer  terms 
as  a  function  of  time.  It  appears  that  the  subharmonic  mode  extracts  most  of 
its  energy  from  the  mean  flow  and  that  there  is  little  net  energy  transfer 
between  it  and  the  fundamental  mode.  In  addition,  its  average  growth  rate 
differs  little  from  that  in  the  absence  of  the  fundamental,  which  is  a  =  0.14. 
Thus,  the  presence  of  the  saturated  two-dimensional  fundamental  does  not  turn 
off  the  subharmonic  mode,  and  the  growth  rate  of  this  latter  mode  is  close  to 
that  of  the  fundamental  two-dimensional  instability.  These  results  imply  that 
even  a  small  subharmonic  perturbation  will  quickly  achieve  finite  amplitude 
after  the  fundamental  mode  saturates,  unless  the  amplitude  of  the  fundamental 
is  artificially  amplified  by  forcing.  In  this  simulation,  the  subharmonic  mode 
saturates  at  t  =  90  at  which  time  the  growth  rate  becomes  negative.  It 
should  be  noted  here  that  in  attempting  to  compare  these  results  for  growth 
rates  of  the  fundamental  and  subharmonic  modes  with  experiments,  it  is  neces¬ 
sary  to  account  for  the  dispersion  of  the  subharmonic  modes  (cf.,  for  example, 
figure  2b  in  Ko  4  Huerre  1984),  which  is  present  in  the  spatially  growing  but 
not  in  the  temporally  growing  mixing  layer.  Also,  in  comparing  these  growth 
rates  with  those  predicted  by  linear  theory,  the  growth  of  momentum  thickness 
of  the  mixing  layer  over  the  course  of  the  simulation  should  be  taken  into 


cons iderat ion . 
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Corcos  &  Sherman  (1934)  find  that  the  presence  of  the  fundamental  inhibits 
the  subharmonic  growth,  while  Pierrehumbert  &  Widnall  (1982)  find  an  enhance¬ 
ment  of  growth.  We  find  that  the  growth  rate  of  the  subharmonic  is  modulated 
with  a  period  related  to  the  oscillation  time  scale  of  the  nonlinear,  quasi- 
equilibrium  fundamental  mode.  However,  the  net  effect  on  the  growth  of  the 
subharmonic  due  to  the  fundamental  is  to  decrease  q  slightly,  from 

about  0.16  to  0.14.  While  these  conclusions  do  agree  with  those  obtained  by 
Kelly  (1967)  using  perturbation  theory,  they  show  that  the  effect  is  quite 
small. 


4.  Three-dimensional  instabilities 

The  saturated  two-dimensional  flow  is  also  subject  to  three-dimensional 

instabilities.  While  the  laminar  mean  flow  is  inviscidly  unstable  only  for 

2  2 

cx  +  3  <1,  the  finite-amplitude  two-dimensional  flow  can  be  unstable  for 

large  G  at  high  Reynolds  numbers.  In  figure  8,  we  plot  the  average  three- 
dimensional  growth  rate  versus  3  of  the  three-dimensional  disturbance 
energy, 

e3d  '  £  •  (;- 

m 

for  various  Reynolds  numbers  when  a  *  0.4  for  a  three-dimensional  linear 
perturbation  to  a  saturated,  two-dimensional  qu3si-equi librium  flow.  Analysis 
of  these  results  suggests  the  conjectures  that,  as  R  increases,  for  a  fixed  3, 
CJD  aPProdC^es  a  finite  limit  ;bo  the  secondary  instability  discussed  is 
inviscid  in  character)  and  that  the  instability  turns  off  for 
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Although  the  mean  flow  tanh  z  is  both  viscously  and  inviscidly  stable  for 
(1  >  1 ,  the  saturated,  two-dimensional  disturbed  flow  is  strongly  unstable  at 
these  scales,  with  disturbances  growing  at  rates  near  those  of  the  inviscid 
two-dimensional  fundamental  instability,  as  shown  in  figure  9.  When  the  two- 
dimensional  modes  saturate,  the  three-dimensional  modes  can  achieve  finite 
amplitudes  on  convective  time  scales  and  thereby  modify  significantly  the  later 
evolution  of  the  flow. 

Like  the  two-dimensional  modes,  the  three-dimensional  modes  evolve  into 
vortical  states  at  large  amplitudes.  Typically,  these  are  counterrotating, 
strearawise  vortices,  or  "ribs",  that  develop  in  the  braids  between  the  two- 
dimensional  vortex  cores.  In  figure  10,  we  show  a  three-dimensional  perspective 
plot  of  surfaces  having  a  value  of  503  of  the  peak  of  the  sum  of  the  absolute 
values  of  all  three  vortieity  components.  The  data  are  from  a  simulation  with 
A1q  -  0.22,  A^2  q  ■  0,  ■  0.003,  R  -  28,  a  »  0.4446,  and  3  »  0.8892  at 

t  *  12.  The  saturated,  two-dimensional  vortex  cores  are  shown  spanning  the 
domain,  and  the  two  pairs  of  counterrotating,  strearawise  vortices  (3  =  2a), 
which  have  evolved  from  a  low-amplitude,  linear  perturbation,  have  developed  in 
the  high-strain  braid  region.  These  are  very  similar  to  the  structures  seen 
experimentally  by  Breidenthal  (1981)  and  Bernal  (1981)  and  modeled  by  Lin  4 
Corcos  (1984);  they  will  be  discussed  in  more  detail  in  §5.  In  this  and  subse¬ 
quent  simulations,  the  choice  $  *  tt/2  'equation  (3.1)]  was  made.  The  choice 
<J>  •  Nil  appeared  to  produce  anomalous  behavior  in  the  evolution  of  the  ribs, 
somewhat  analogous  to  the  two-d fcnens iona  1  shredding  interaction  when  9  *  Nrr. 

The  cutoff  as  3  ■*  3  (4.2)  is  ma.nly  due  to  increasing  viscous  damping, 

cri  t 

as  nonlinear  transfers  vary  less  rapidly  with  3.  This  point  is  suggested  by 
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Here  Yw  involves  the  nonlinear  interaction  of  the  three-dimensional  and  mean 
M 

flows,  and  Y->D  the  interaction  of  the  three-dimensional  and  two-dimensional 
energy.  The  jitter  in  the  plots  is  due  to  numerical  inaccuracy  in  the  evalua¬ 
tion  of  the  nonlinear  transfer  terms,  so  only  the  trends  in  the  data  should  be 
considered  significant.  It  is  suggested  from  the  data  in  figure  lia  that,  for 

Q  <  Q  .  ,  Y  and  <<  Y..  asymptotically  so  that  the  three-dimensional  mode 
crit*  'v  2D  M 

derives  its  energy  from  the  mean  flow  with  the  two-dimensional  disturbance 
acting  as  a  catalyst  for  this  transfer. 

On  the  other  hand,  the  results  plotted  in  figure  lib  show  that  when 
B  ~  ^cr^ci  Yv  is  quite  significant.  The  three-dimensional  instability  seems  to 
be  turned  off  at  a  large  cross-stream  wavenumber  ft  by  increased  dissipation 
rather  than  by  any  significant  qualitative  change  in  nonlinear  transfers  from 
the  mean  and  two-dimensional  components. 

The  nature  of  the  competition  between  two-dimensional  pairing  and  three- 
dimensional  instability  is  illustrated  by  the  results  plotted  in  figures  12  and 
13.  In  both  figures,  we  have  plotted  the  results  of  runs  with  R  *  400,  a®  0.4 
and  (1  *  0.2.:  Figure  12  shows  the  evolution  of  the  instabilities  when  the 
initial  three-dimensional  perturbation  is  much  larger  than  the  subharmonic 
mode.  In  this  case,  the  pairing  instability  is  nearly  unaffected  by  the  three- 
dimensional  instability  before  finite  amplitudes  are  reached.  In  figure  13, 
the  initial  conditions  are  chosen  so  that  the  subhar-or.ic  mode  perturbation  vs 


-im¬ 


pairing  process  and  rollup  inhibit  the  three-dimensional  instability  near 
t  =  72,  where  0^  actually  becomes  negative.  Note  that  the  growth  rate  of 
E-d  is  identical  for  the  two  cases  (figures  12  and  13)  until  ^  reaches 

finite  amplitude. 

5 .  Instability:  Dependence  on  initial  conditions 

The  flows  that  develop  from  the  three-dimensional  secondary  instability  do 
not  appear  to  saturate  in  ordered  states  like  those  of  the  fundamental  and 
subharmonic  two-dimensional  instabilities.  However,  the  presence  of  additions 
pairing  modes  can  significantly  enhance  the  overall  coherence  of  the  flow. 
Experiments  have  shown  that  extremely  low  levels  of  forcing  can  generate 
changes  in  the  flow  of  order  1  (e.g.,  Gutmark  &  Ho  1983).  This  has  raised 
concerns  that  unintended  forms  of  weak  excitation  due  to  natural  resonances  in 
an  experimental  apparatus,  or  pressure  feedback  effects,  could  be  modifying  th 
evolution  cf  the  flow  field  (Hussain  1983b).  An  important  difference  between 
our  numerically  simulated  temporally  growing  flows  and  the  spatially  growing 
flows  in  laboratory  experiments  is  that  influences  of  downstream  events  on  the 
earlier  stages  of  evolution  are  possible  in  the  latter  but  not  the  former. 

In  figure  14a,  we  plot  results  that  show  the  effect  of  the  absence  of  a 

subharmonic  pairing  mode  on  the  evolution  of  the  modal  energies.  In  the 

absence  of  the  subharmonic,  the  fundamental  quickly  rolls  up,  a  process  that 

inhibits  the  growth  of  the  low-amplitude  three-dimensional  disturbance.  Once 

0 

the  fundamental  reaches  its  saturated  quasi-equilibrium  state  (t  -  20),  the 
three-dimensional  modes  resume  their  rapid  growth.  By  about  t  *  50,  the  three 
dimensional  modes  dominate  the  flow  field.  The  perspective  plot  given  in 
figure  14b  shows  this  domination  even  at  later  times.  Further  growth  of  the 
three-d iraens iona 1  perturbation  energy  is  inhibited  by  the  collapse  of  the  mixi 


layer  due  to  the  absence  of  the  subharnon ic .  With  the  subt-  lr-nonic  present,  the 
evolution  of  the  three-dimensional  modes  is  dramatically  different.  The 
results  plotted  in  figure  14c  are  from  a  simulation  identical  to  the  previous 
one  except  for  the  inclusion  of  the  subharmonic  mode.  The  three-dimensional 
modal  growth  (at  amplitudes  in  the  linear  range)  is  now  slowed  both  by  the 
fundamental  rollup  (t  s  5-10)  and  by  the  subharmonic  pairing  (t  =  20-35). 

Thus,  the  flow  is  more  coherent  than  in  the  absence  of  the  subharmonic  (compare 
figures  14b  and  14d).  Once  the  subharmonic  reaches  its  saturated  state 
(t  *  35),  the  three-dimensional  growth  rate  increases  substantially. 

As  was  pointed  out  in  54,  the  two-dimensional  and  three-dimensional 
unstable  modes  grow  at  very  similar  rates  when  all  disturbance  amplitudes  are 
in  the  linear  range.  In  figure  15,  we  plot  the  results  of  a  simulation  in 
which  both  fundamental  and  three-dimensional  modes  were  introduced  at  approxi¬ 
mately  equal  amplitudes,  well  down  into  the  linear  range.  In  this  case,  the 
three-dimensional  mode  disrupts  the  rollup  and  saturation  of  the  fundamental 
mode,  so  that  its  peak  amplitude  is  an  order  of  magnitude  smaller  than  without 
the  three-dimensional  mode  present.  As  shown  in  figure  15b,  the  presence  of 
the  large  three-dimensional  instability  substantially  reduces  the  spanwise 
coherence . 

To  determine  whether  the  model  used  to  initiate  the  three-dimensional 
instabilities ■ in  the  simulations  described  thus  far  was  realistic,  we  performed 
several  simulations  initialized  with  a  broad  range  of  three-dimensional  modes. 
We  used  an  uncorrelated,  random-phase  velocity  field  having  a  Gaussian-shaped 
energy  spectrum.  This  field  was  convolved  with  a  function  so  that  the  relative 
turbulence  intensity  levels  were  consistent  with  those  of  experimental  mixing 
layer  data.  However,  the  initial  peak  intensity  level  was  about  3  orders  of 


magnitude  below  the  experimental  values,  so  the  initial  disturbance  growth  was 
in  the  linear  regime. 

In  figure  16a,  we  plot  the  evolution  of  the  flow  field  with  broad-band 
initial  excitation.  Neither  fundamental  nor  subharmonic  two-dimensional  linear 
eigenfunction  modes  were  explicitly  included  in  the  initial  conditions,  although 
there  was  energy  in  the  corresponding  wavenumbers  defined  by  the  random  ini¬ 
tialization  process.  Nonetheless,  E^q  and  q  Srow  ve*T  rapidly  initially, 

with  a1Q  *  0.11  and  q  *  0.13  at  t  *  10  (compared  with  q  *  0.19  and 

*  0.14  for  linear  inviscid  modes).  E^,  which  is  the  energy  in  all  velocity 
components  not  having  k^  *  0,  has  a  growth  rate  of  *  0.12  at  t  *  10.  When 
the  modal  components  of  Ej  reach  finite  amplitude  (at  t  -  35),  cr^  decreases 
sharply  and  the  flow  field  is  characterized  by  a  chaotic,  turbulent-like  velo¬ 
city  field  in  the  center  of  the  mixing  layer  (figure  16b).  The  flow  remains  in 
this  chaotic  state  until  the  subharmonic  reaches  nonlinear  amplitude  and  begins 
to  roll  up.  As  the  subharmonic  reaches  its  saturation  amplitude  (t  *  90), 
the  nature  of  the  three-dimensional  velocity  field  evolves  from  a  highly  chaotic 
state  to  one  with  coherent,  large-scale  structures  (compare  figures  16b  and 
16c),  although  remains  approximately  constant.  Thus,  there  is  a  complex 
interaction  between  the  two-  and  three-dimensional  instabilities.  The  coherent, 
two-dimensional  modes  significantly  enhance  Che  growth  of  the  mixing  layer 
thickness,  creating  a  larger  region  for  the  ultimate  expansion  of  three- 
dimensional  instabilities.  During  their  rollup  and  pairing  processes,  however, 
the  coherent,  two-dimensional  modes  tend  to  reduce  significantly  the  growth 
rates  of  the  low-amplitude  three-dimensional  modes. 

The  inhibition  of  the  three-dimensional  modal  growth  rate  by  the  rollup  and 
pairing  of  the  two-dimensional  modes  is  a  function  of  both  the  amplitude  and 


tl'.e  wavenumber  of  Che  three-dimensional  modes.  For  low-amplitude  three- 
dimensional  modes,  in  the  linear  range,  coherent  rollup  or  pairing  can  actually 
be  stabilizing.  In  figure  17,  we  plot  the  evolution  of  the  modal  energies  as 
functions  of  time  for  a  simulation  with  the  same  initial  conditions  as  in 
figure  16  but  with  q  “  0*0014.  As  seen  in  this  figure,  when  the  funda¬ 

mental  and  subharmonic  reach  nonlinear  amplitudes  (t  =  25),  the  three- 
dimensional  growth  is  stopped  completely.  At  higher  three-dimensional  modal 
amplitudes,  however,  the  rollup  may  have  little  or  no  effect  on  For 

example,  in  figure  16a,  cr^  changes  very  little  in  the  presence  of  the 
saturating  subharmonic  between  t  *  60  and  100.  During  the  three-dimensional 
stabilization  shown  in  figure  17,  the  low  3  modes  continue  to  grow  while  the 
high  fl  modes  decay.  This  is  not  unreasonable,  since  as  the  scale  of  the  mixing 
layer  grows  due  to  the  two-dimensional  pairing,  the  relative  scale  of  the 
spanwise  instability  also  changes.  As  shown  in  figure  8  for  three-dimensional 
perturbations  to  two-dimensional  saturated  modes,  depends  strongly  on 

(J.  A  similar  effect  is  seen  for  E„  in  the  nonlinear  range  in  figure  16a. 

This  would  perhaps  be  manifest  in  the  formation  of  larger  ribs  after  the  pair¬ 
ing  of  the  two-dimensional  modes.  One  mechanism  that  may  have  a  significant 
influence  on  the  suppression  of  the  three-dimensional  modal  growth  is  the 
temporal  variation  of  the  strain  field  in  the  braids  between  the  coherent, 
two-dimensional  vortices.  In  the  early  stages  of  rollup,  a  very  high  strain 
develops  in  the  braids.  This  has  a  tendency  to  stretch  the  ribs,  intensifying 
the  streamwise  vorticity.  As  the  two-dimensional  modes  approach  saturation, 
however,  the  strain  race  decreases  substantially,  so  that  this  vortex  stretch¬ 


ing  mechanism  is  weakened. 
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o •  The  secnndar v  instability  in  a  turbulent  flow 

We  have  performed  high-resolution  (64x64x64  mode)  simulations  of  a  fully 
turbulent  mixing  layer,  and  it  is  instructive  to  relate  the  evolution  of  these 
flows  to  Che  class  of  instabilities  discussed  so  far.  These  simulations  were 
performed  on  a  computational  domain  sufficiently  large  to  allow  two  complete 
pairings  (periodicity  length  8u/a).  The  initialization  procedure  was 
similar  to  that  discussed  in  §5,  but  the  amplitudes  of  the  initial  fields  were 
higher.  Details  of  the  numerics  are  given  in  Riley,  Metcalfe  &  Orszag  (1985). 

The  spanwise  vorticity  field  after  two  complete  pairings  shows  clear 
evidence  of  large-scale  structures  (see  figures  18a  and  b).  A  comparison  of 
the  vorticity  plots  at  two  different  spanwise  locations  indicates  a  strong 
spanwise  coherence  for  this  particular  realization,  although  details  of  the 
structures  are  different.  As  previously  noted,  the  secondary  instabilities  in 
the  mixing  layer  flow  are  characterized  by  streamwise,  counterrotating  vortices 
that  tend  to  form  in  the  braids.  Figure  19  is  a  contour  plot  of  in  a 
plane  at  the  middle  of  the  mixing  layer  (z  *  0).  The  solid  and  dashed  lines 
indicate  positive  and  negative  vorticity,  respectively.  Figure  20  is  a  similar 
plot  at  a  streamwise  location  B/4  from  the  left  boundary  in  figure  19.  These 
two  plots  clearly  indicate  the  presence  of  such  vortices,  although  they  are 
irregularly  spaced. 

One  of  the  best  laboratory  visualizations  of  coherent,  three-dimensional 
structures  in  a  turbulent  mixing  layer  was  performed  by  3ernal  &  Roshko  (cf. 
Bernal  1981).  Using  laser  fluorescein  dye  techniques,  they  were  able  to 
illuminate  the  flow  at  a  fixed  streamwise  location.  The  most  striking 
characteristic  of  these  photographs  is  the  appearance  of  mushroom-shaped 
features  on  the  braids  between  the  two-dimensional  vortex  cores  (figure  21a). 

We  have  been  able  to  simulate  this  technique  by  employing  a  numerical  cod'1 


developed  Co  study  a  chemically  reacting  mixirg  layer  (Riley  e t  a  1 .  19.-15)  in 
which  the  advection  diffusion  equations  for  a  binary  chemical  reaction  are  solved 
along  with  the  Navier-Stokes  equations.  Figure  21b  is  a  contour  plot  of  the 
concentration  of  one  of  the  chemical  species  on  the  same  plane  as  in  figure  20, 
and  figure  21c  corresponds  to  the  concentration  of  the  other  chemical  species. 

A  comparison  of  figures  20  and  21b  shows  that  the  counterrotating  vortex  pairs 

tend  to  pump  fluid  through  the  braid  between  their  cores,  increasing  the  reac-  , 

tion  surface  area  and  creating  the  mushroom  shaped  structures  in  the  flame 

front,  which  is  defined  by  the  region  of  overlap  between  figures  21b  and  21c. 

Such  features  were  also  noted  in  the  model  proposed  by  Lin  &  Corcos  (1980. 

The  structure  of  the  streamwise,  counterrotating  vortex  tubes  is  made 
clearer  by  isolating  the  streamwise  vorticity  component  of  the  flow.  Figure  22 
is  a  three-dimensional  perspective  plot  of  surfaces  at  which  equals 

50%  of  its  peak  value.  This  figure  was  from  the  same  realization  and  at  the 
same  time  as  figures  18-21.  The  large-scale,  spanwise-coherent  structures  do 
not  show  up  in  this  plot  since  they  consist  mainly  of  spanwise  vorticity, 

Wy.  Comparison  with  figures  18a  and  18b  shows  that  the  ribs  do  form  on  the 
braid  between  the  large-scale  two-dimensional  vortex  cores.  This  structure  is 
consistent  with  the  model  proposed  by  3ernal  (1981)  although  the  irregular 
spacing  of  the  ribs  suggests  that  the  modification  of  this  model  proposed  by 
Hussain  (1933a)  is  more  realistic.  The  effect  of  increased  coherence  of  the 
t wo-dimens ional  pairing  modes  on  the  rib  structure  is  shown  in  figure  23,  which 
is  from  a  simulation  like  the  previous  one  but  to  which  two-dimensional  modes 

have  been  added  in  the  initial  conditions:  A, „  *  0.1,  A. a  0.06, 

i  J  4.  /  Z  f  J 

A.  nm  0.025  [equation  (3.1)].  The  resulting  ribs  are  more  coherent  and 
i  /  •* ,  u 

more  aligned  in  the  streamwise  direction.  As  was  the  case  with  the  simulation 
in  figure  lb,  the  presence  of  the  pairing  mode  tends  to  increase  the  coherence 
of  the  three -dimensional  perturbation  field. 
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The  extreme  sensitivity  of  the  ribs  to  the  initial  or  upstream  flow 
conditions  makes  direct  quantitative  comparisons  between  the  simulations  and 
laboratory  experiments  difficult.  In  the  simulations  we  have  performed  so  far, 
there  have  been  significant  variations  in  the  amplitudes  and  spanwise  spacing 
of  the  ribs.  Likewise,  in  the  experiments  of  Bernal  (1981),  there  was 
substantial  scatter  in  the  measurements  of  the  rib  spacing.  In  addition,  he 
found  that  the  spanwise  position  of  the  ribs  appeared  to  be  related  to 
disturbances  originating  upstream  in  the  settling  chamber.  Our  stability 
analysis  (figure  8)  has  shown  that  there  is  a  broad  range  of  spanwise 
wavenumbers  that  are  unstable.  Some  representative  values  for  the  simulation 
shown  in  figures  18-22  are  as  follows:  the  peak  spanwise  vorticity  normalized 
by  the  peak  mean  velocity  gradient  is  about  2,  while  the  peak  strearawise 
vorticity  is  slightly  higher,  about  3.  The  spanwise  vorticity  amplitudes  are 
consistent  with  filtered  experimental  data  (Metcalfe,  Hussain  &  Menon  1985), 
while  the  streamwise  vorticity  amplitudes  are  somewhat  higher  than  i.n  other 
runs.  The  rib  spacing  (estimated  from  figure  19)  is  about  the  same  as  the 
wavelength  of  the  most  unstable  fundamental  two-dimensional  mode.  This  is  in 
the  range  of  values  reported  by  Bernal  (1981).  A  more  detailed  analysis  of 
these  simulations,  using  experimental  data  to  refine  these  comparisons,  is  now 
in  progress. 

7 .  Discussion 

We  have  shown  that  small-scdle  three-dimensional  instabilities  like  those 
previously  studied  by  Orszag  &  Patera  (  1980,  1983)  and  P i er rehumber t  &  Widnall 
(1982)  exist  in  free  shear  flows  and  that  these  instabilities  persist  to 
moderately  low  Reynolds  numbers.  It  is  now  clear  that  these  modes  can  be 
responsible  for  the  initial  development  of  three- d  -.men-. .  ->na!  itv  in  these  she  l: 
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flows.  The  dynamics  of  these  three-d  i:v*ns  ion.,  1  instabilities  is  sir.il.ir  t  a 
that  of  the  t hree-d imens iona 1  i ns t ab i 1 1 1 ics  i n  wall-bounded  shear  flows.  At 
high  amplitudes,  these  instabilities  manifest  themselves  mainly  as  counter¬ 
rotating,  streamwise  vortices,  or  ribs,  that  form  on  the  braids  between  the 
spanwise-coherent ,  two-dimensional  pairing  modes;  they  are  responsible  for  the 
generation  of  the  mushroom-shaped  features  seen  in  laboratory  experimental 
visualizations  (Bernal  1981).  While  the  instabilities  share  some  features  of  a 
classical  inflectional  instability,  including  phase  locking  with  the  funda¬ 
mental  vortex,  inflectional  instabilities  are  preferentially  two-dimensional, 
whereas  the  present  instabilities  are  not. 

The  rollup  and  pairing  of  the  two-dimensional  modes  has  a  stabilizing 
effect  on  the  higher  spanwise  modes  and  on  the  overall  three-dimensional  growth 
rate  when  the  amplitude  of  the  three-dimensional  modes  is  small,  while  the 
absence  of  pairing  (saturation)  can  enhance  the  three-dimensional  growth  rate. 
However,  in  the  absence  of  pairing  or  rollup,  the  three-dimensional  instabil¬ 
ities  can  reach  a  chaotic,  saturated  state  from  which  significant  further 
growth  is  not  possible  without  additional  pairing.  The  suppression  of  the 


low-amplitude  three-dimensional  instabilities  by  pairing  could  explain  the 
strong  two-dimensionality  of  the  flow  near  the  splitter  plate  in  many  labora¬ 
tory  experiments.  Once  the  three-dimensional  modes  reach  a  finite  amplitude 
and/or  the  Reynolds  number  increases  with  downstream  distance,  the  growth 
suppression  effect  is  reduced,  and  the  flow  becomes  mere  three-dimensional . 
Mixing  can  be  enhanced  initially  by  coherent  forcing  of  the  mixing  layer  to 
saturation  (Hussain  S  Zaraan  1978  ;  Oster  &  Wyg-anski  1 3  3  2 ;  Ho  &  Huang  19S2), 
which  will  enhance  the  growth  of  the  three-dimensional  modes  to  a  chaotic 
state.  But  this  is  done  at  the  expense  of  reducing  the  growth  of  the  mixing 


l.iver  momentum  thickness,  which  is  due  primarily  to  the  pairing  nrorejs.  Th  ;  , 


the  maximization  o  t  pr'duct  formation  in  j  rrart  irn;  mixing  layer  will  ri'rrur- 
balancing  the  increases  in  the  flame  front  area  generated  by  the  chaotic 
three-dimensional  flow  with  the  increases  Jue  to  the  more  rapid  mixing  layer 
growth  caused  by  the  presence  of  additional  pairing  modes. 

The  results  of  these  simulations  suggest  that,  with  respect  to  the  growth 
of  three-dimensional  d isturbances ,  there  are  several  important  flow  states 
possible  in  an  evolving  mixing  layer.  First,  there  is  pairing  and  rollup  of 
the  two-dimensional  modes,  which  is  characterized  by  a  suppression  of  low- 
amplitude  three-dimensional  modal  growth  at  low  Reynolds  numbers.  Second, 
there  is  the  saturated,  two-dimensional,  quasi-equilibrium,  non-pairing  state 
which  is  highly  unstable  to  three-dimensional  perturbations.  Finally,  there 
are  chaotic,  three-dimensional  states  characterized  by  a  lack  of  spanwise 
coherence  from  which  only  moderate  three-dimensional  growth  can  occur.  It  is 
from  these  states  that  more  rapidly  growing,  large-scale,  two-dimensional  mod 
can  eventually  emerge  and  reorganize  the  flow  in  a  manner  consistent  with  tha 
suggested  by  Staquet  &  Lesieur  (1986). 

It  seems  that  the  mechanics  of  transition  in  the  free  shear  flows  studied 
here  may,  in  a  sense,  be  rather  more  complicated  than  in  the  case  of  wall- 
bounded  shear  flows.  In  the  latter  case,  linear  instabilities  are  often 
viscously  driven  and,  therefore,  weak,  so  they  cannot  be  directly  responsible 
for  the  rapid  distortions  characteristic  of  transition.  On  the  otner  hand, 
free  shear  flows  are  subject  to  a  variety  of  inviscid  instabilities,  so  there 
may  be  many  paths  to  turbulence!  We  have  shown  that  the  choice  of  paths  in  a 
individual  flow  may  depend  an  the  results  of  cc-rnetiticn  between  fundament  -•  I  , 
subharmonic,  and  three-dimensional  instabilities,  all  of  which  are  convict  ;v>- 
driven  and,  therefore,  strong  with  comparable  gr'wrh  rote?.  Thus,  the  •»-.•- 1 - 
t  ion  of  free  shear  flows  in  transitional  regime  depend  ;  s  i  gni  l  ie  .1:.:  1  v  .  n 
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pea t  history  of  the  flow.  includins  it*  relative  onpillnde:.  of  .11  ct.pt 
nodes ,  the  aechaniso  of  their  d«neratlo».  an  1  toe  external  environment  1 

the  flow  Is  embedded. 
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Jrowfh  rates  (In  c)  of  tbs  0rr-.3omner  f<*  1.!  eigenfunctions 
£^r  the  nixing  layer  n  q  ( z. )  *  Ug  tanh  (r,'tj)t 


x-vavenumber  a 


Number  of  Chebyshev  Polynomials  (P+1. 


Hyperbolic  Map  (2.5, 


1.534 

0.959 

0.635 

0.612 

0.539 

0.202 


1.375 

0.320 

0.614 

0.598 

0.597 

0.526 


1.238 

0.746 

0.605 

0.597 

0.597 


0.579 

0.383 

0.344 

0.324 

0.115 


0.596  Stt 


0.501 

0.360 

0.342 

0.342 

0.322 
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Algebraic  Map 


0.457 

0.351 

0.342 

0.342 

0.342 

0.321 


0.160 

0.141 

0.137 

0.041 

S 

s 


0.5 

0.699 

0.588 

0.599 

0.345 

0.346 

0.342 

0.131 

0.133 

0.137 

1 

0.591 

0.599 

0.597 

0.344 

0.342 

0.342 

0.137 

0.137 

0.137 

2 

0.600 

0.597 

0.597 

0.342 

0.342 

0.342 

0.136 

0.137 

0.137 

4 

0.597 

0.597 

0.597 

0.325 

0.342 

0.342 

0.371 

0.136 

0.137 

8 

0.542 

0.597 

0.597 

0.009 

0.322 

0.342 

S 

0.043 

0.136 

T  Here  the  Reynolds  number  is  Uq6^/v  ■  100  and  the  eigenvalue  is  the  complex 
wave  speed  c  for  a  temporal  mode  of  the  form  t(z)e^^x'ct^.  For  the  most 
rapidly  growing  code  listed  here,  Re(c)  =  0. 


ft  S  indicates  that  all  modes  are  stable  with  the  indicated  parameter  values 
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Figure  Copt  ions 

FIGURE  1.  A  plot  of  E^q(c)  versus  t  for  runs  with  A,  ^  *  A,  ^  3  0  and 

A^q  ■  0.5,  0.25,  0.125,  0.01.  Here  the  Reynolds  number  is 

R  *  400,  the  spectral  cutoffs  in  (2.4)  are  M  *  8 ,  N  ■  1 ,  P  *  32 
(resolution  8x1x32  with  no  subharraonic  modes),  the  x-wavenumber  is 
a  »  0.4,  and  the  time  step  is  At  *  0.02.  Note  that  the  flow 
saturates  into  a  vortical  state  nearly  independent  of  the  initial 
perturbation.  Before  such  saturation  occurs,  the  perturbation 
grows  linearly  like  an  Orr-Sommerfeld  eigenfunction. 

FIGURE  2.  Contour  plots  of  spanwise  (y)  vorticity  for  the  mixing  layer  at 

t  ■  8  and  t  *  32  with  R  a  83,  A^,  Q  *  0,  A1Q  *  0.2, 

M  •  P  •  64,  and  N  *  1. 

FIGURE  3.  Plots  of  the  evolution  of  E^g(t)  and  the  two-dimensional 

subharmonic  mode  energy  E g(t).  Here  R  *  400,  A^g  *  0.25, 

A 1 / 2  g  *  3x10  4,  M  •  8,  N  »  l,  P  ■  32,  a  *  0.5,  At  -  0.02, 
and  8  *  7r/2a  (phase  difference  between  the  two  modes). 

FIGURE  4.  Spanwise  vorticity  contours  at  t  “  0,  8,  16,  and  24,  during  a 
vortex  pairing  run  with  R  »  83,  A^g  »  0.20,  q  “  0.14, 

M  ■  64,  N  *  1,  P  ■  64  (resolution  64x1x64),  a  *  0.4446, 

At  ■  0.05,  and  9  *  7r/2a. 

FIGURE  5.  Velocity  vector  plots  of  the  vortex  pairing  run  shown  in  figure  4. 

FIGURE  6.  (a  -  d)  Plots  of  the  Reynolds  stress  at  different  times  for  the 
runs  with  both  fundamental  and  subharmonic  present  (figure  4). 

(e)  A  plot  of  the  Reynolds  stress  with  and  without  pairing  at 
different  times. 
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FIGURE  7.  A  plot  of  the  components  y  ,  Y,n,  Y  [see  (3.3)!  of 

M  *-  J  V 

the  growth  rate  j7  Q  subharmonic  mode  amplitude  as 

functions  of  time  for  a  run  with  R  =  200  ,  A^q  3  0.25, 

A j q  -  3x10  ^ ,  M  3  16,  H  3  1,  P  3  32,  a  3  0.43,  and 
At  ■  0.01. 


FIGURE  8.  A  plot  of  the  computed  three-dimensional  growth  rate  [see 

(4.3)]  as  a  function  of  the  spanwise  wavenumber  II  for  various 
Reynolds  numbers.  Here  a  3  0.4. 


FIGURE  9.  A  plot  of  the  computed  subharmonic  growth  rate  0^/2  q  and 
three-dimensional  growth  rate  G^  as  a  function  of  a  at 
R  *  400  and  13  *  0.3. 


FIGURE  10.  A  three-dimensional  perspective  plot  of  surfaces  having  a  value 

equal  to  50%  of  the  peak  of  the  sum  of  the  absolute  values  of  all 
three  vorticity  components  for  a  run  with  R  3  56  (at  t  *  12), 


M  3  N 


64,  a  -  0.4446,  At  *  0.05,  A 


0.22,  A 


Au  -  0.003,  .3  *  0.8892  at  t  -  12. 
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FIGURE  11.  A  plot  of  the  components  y  ,  y y  [see  (4.3)]  of  the 

4*  *m  U  V 

three-dimensional  growth  rate  G^  as  functions  of  time  for 
R  -  400,  a  -  0.4,  A1q  »  0.25,  An  -  10~\ 

(a)  3  3  4. 

(b)  3-6. 


FIGURE  12.  A  plot  of  the  evolution  of  the  energies  E.  _ ,  E.  ,,,  E, 

1U  1  /  *. ,  U  3  u 

versus  t  for  a  run  with  R  =  400,  a  3  0.4,  3  3  0.2,  M  3  8,  N  3  4, 
P  3  32,  and  initial  conditions  A^q  3  0.25,  A^^  q  *  3x10  , 

A^  3  10  ^ .  The  thr^e-ciimens iona  1  mode  initially  dominates  the 
subharmonic  mode. 
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12,  except  th.it  the  initial  conti  it  ions  ire 


°,:5,  A 1 / 2 , 0 


4x10  \  A ,  ^  3  3.3x10  h  Here, 


from  figure  12.  Note  that  the  linear  growth  oi  the  three-d imens iona  1 
mode  is  stopped  when  the  subharmonic  reaches  a  nonlinear  amplitude. 

FIGURE  14.  (a)  A  plot  of  the  evolution  of  the  modal  energies  as  functions  of 
time.  Eg  3  Eq^  ♦  ,  i.e.,  the  three  lowest  spanwise 

nodes  with  the  same  x-vavelength  as  the  fundamental.  R  3  100, 
a  -  0.4446,  A^q  3  0.22,  A^2  Q  *  0 ,  and  A^  3  10-5. 

(b)  A  three-dimensional  perspective  plot  of  the  vorticity  field  as 
in  figure  10  at  t  *  96. 

(c)  Same  as  figure  14a  but  \  •  0.14. 
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(d)  Same  as  figure  14b  but  A^2  q  3  0.14. 


FIGURE  15.  (a)  Same  as  figure  14a  but  A^q  ■  10  . 

(b)  Same  as  figure  14b  but  A^g  e  10  ^  at  t  **  48. 


FIGURE  16.  (a)  A  plot  of  the  evolution  of  the  modal  energies  as  functions  of 
tine.  Er  *  energy  in  all  modes  with  k  ^0.  R  *  100.  Random 

-  y 

noise  initial  field  with  peak  rms  velocity  =  0.01. 

(b)  A  three-dimensional  perspective  plot  of  the  vorticity  field  as 
in  figure  10  at  t  *  48. 

(c)  Same  as  figure  16b  but  at  t  3  96. 

FIGURE  17.  A  plot  of  the  evolution  of  the  modal  energies  as  functions  of  time. 
Er  3  energy  in  all  modes  with  k^  f  0.  R  3  100, 
a  -  0 . 44h6  ,  A10  =  0,  A1/2|0  3  0.0014,  and  Au  3  10  . 


FIGURE  18.  Spanwise  vorticity  at  y  *  4r/a  (figure  18a)  and  6ir/a  (figure  13b). 
The  computational  doc^in'size  B  is  8 tt/  a.  R  *  154  and  t  3  72. 
Random  noise  initial  field  with  peak  rms  velocity  -  0.13.  The 
initial  Reynolds  number  is  R  3  28. 
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FIGURE  19.  Strea-.nwise  vor'  icity  at  z  =  )  tor  the  sa-e  run  as  in  figure  1.3. 

FIGURE  20.  Screamwi.se  vorticity  at  x  a  2T/a  for  the  same  run  as  in 

figure  IS. 

FIGURE  21.  (a)  Laser  sheet/fluorescein  dye  visualization  of  Che  braid  of  a 
turbulent  mixing  layer  in  water  (from  Bernal  1981). 

(b)  Contour  plot  of  the  species  concentration  field  at  the  same 
time  and  location  as  in  figure  20. 

(c)  Contour  plot  of  the  second  species  concentration  field  at  the 
same  time  and  location  as  in  figure  21b. 

FIGURE  22.  Three-dimensional  perspective  plot  of  surfaces  at  which  Iuj^I 

equals  50%  of  its  peak  value  for  the  same  run  as  in  figure  18. 

FIGURE  23.  Plot  of  iu  I  as  in  figure  22  but  for  a  simulation  to  which 

two-dimensional  modes  have  been  added  in  the  initial  conditions: 
A, -  *  0.1,  A.  -  *  0.06,  and  A,  ,.  A  a  0.025  [equation  (3.1)]  . 


PEAK  CONTOUR  VALUE  IS -1.12.  CONTOUR  INTERVAL  IS  0.07.  t-8. 
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